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Abstract. We study the Hartree-Fock model for pseudorelativistic atoms, 
that is, atoms where the kinetic energy of the electrons is given by the pseudo- 
relativistic operator ^(|p|c)2 + (mc^)^ — rric^ . We prove the existence of a 
Hartree-Fock minimizer, and prove regularity away from the nucleus and point- 
wise exponential decay of the corresponding orbitals. 



1. Introduction and results 

We consider a model for an atom with N electrons and nuclear charge Z , where 
the kinetic energy of the electrons is described by the expression -y/dplc)^ -f {mc^Y— 
mc^ . This model takes into account some (kinematic) relativistic effects; in units 
where h = e = m = \, the Hamiltonian becomes 

^ Z 



l<i<j<N 



N 

j=l l<i<j<N 



with T(p) = £;(p) - = \/\pF+a^ - and V{x) = Za/\x\. Here, a is 
Sommerfeld's fine structure constant; physically, a ~ 1/137.036. 

The operator H acts on a dense subspace of the A^-particle Hilbert space Hp = 
Afl]^L^(]R^; C) of antisymmetric functions, where q is the number of spin states. 
It is bounded from below on this subspace (more details below). 

The (quantum) ground state energy is the infimum of the spectrum of H consid- 
ered as an operator acting on Tip: 

E^^{N,Z,a) :=infcr^^(i7) =inf {q(*, «')!*£ Q(i/),(*,*) = 1}, 

where q is the quadratic form defined by H, and Q the corresponding form domain 
(see below); ( , ) is the scalar product in T^f C (M^w • c<?" ) , 

In the Hartree-Fock approximation, instead of minimizing the functional q in the 
entire iV-particle space Tip , one restricts to wavefunctions ^! which are pure wedge 
products, also called Slater determinants: 

1 

with {ui}fLi orthonormal in L'^{M.^;C^) (called orbitals). Notice that this way, 
^ GHp and ||«'||i2(K3«.c,«) = 1. 



4'(xi,cri;x2,(T2; . . . ;xAr,crjv) = det(Mi(xj,crj))j^j^i , (2) 
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The Hartree-Fock ground state energy is the infimum of the quadratic form q 
defined by H over such Slater determinants: 

E^^(N, Z, a) -.^ inf{ q(*, *) | * Slater determinant} . (3) 

For the non-relativistic Hamiltonian, 

N 

H.(iV,Z) = i:{ A, }+ Y. (4) 

j=l ' ' l<i<j<N ' ^ ' 

the mathematical theory of this approximation has been much studied, the ground- 
breaking work being that of Lieb and Simon [TH] ; see also [H] for work on excited 
states. For a comprehensive discussion of Hartree-Fock (and other) approximations 
in quantum chemistry, and an extensive literature list, we refer to [10'. 

The aim of the present paper is to study the Hartree-Fock approximation for the 
pseudorelativistic operator H in 

We turn to the precise description of the problem. The one-particle operator ho = 
T(-iV) - y(x) is bounded from below (by a'^ [(1 - (7rZa/2)2)i/2 - 1]) if and only if 
Za < 2/7r (see [7], [9j 5.33 p. 307], and [25,; we shall have nothing further to say on 
the critical case Za = 2/7r). More precisely, if Za < 1/2, then F is a small operator 
pertubation of T. In fact [7, Theorem 2.1 c)], \\\Ti\-^{T{-iV) + =2. 

As a consequence, ho is selfadjoint with T>{ho) = i/^(M'^; C) when Za < 1/2. It is 
essentially selfadjoint on C^(E3;C9) when Za < 1/2. 

If, on the other hand, 1/2 < Za < 2/tt, then V is only a small form pertubation 
of T: Indeed 5.33 p. 307], 

/ ^^dx<5 / |p||/(p)pdp for /ei/i/2(M3), (5) 

where / denotes the Fourier transform of /. Hence, the quadratic form o given by 

V[u,v] {V^^'^u,V^^^v) for u,v e H^/^{R^;C'') (6) 

(multiplication by 1/^/^ in each component) is well defined (for all values of Za). 
Here, ( , ) denotes the scalar product in L^(M'^;C'^). Let e be the quadratic form 
with domain H'^/'^(R^;C'^) given by 

e[w, v] (^(p)i/2u, E{py/^v) for u,v e //^/^(M^'; C«) . (7) 

By abuse of notation, we write £'(p) for the (strictly positive) operator i?(— iV) = 
V-A + q;-2. Then, using (O and that |p| < E{p), 

V[u,u] < e[u,u] for u e H^^^ (R^iC^) if Za <2/Tr. (8) 

Hence, by the KLMN theorem Theorem X.17], there exists a unique self-adjoint 
operator ho whose quadratic form domain is i/^/^(IR'^; C) such that (with t = 
e-a-i) 

iu,hov) = t[u,v]-v[u,v] for u,v e H^^^iR^iC^) , (9) 

and ho is bounded below by —a~^. Moreover, if Za < 2/7r then the spectrum of 
ho is discrete in [—a~^, 0) and absolutely continuous in [0, cxd) [71 Theorems 2.2 and 
2.3]. 
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As for the A^-particle operator in when Za < 2/7r, ^ implies that the 
quadratic form 

N 

N 

+ J2 (|x,-X,rl/2vI/,|x,-X,|-l/2$) , xP,$e /\ifl/2(R3.c9)^ 

l<i<j<N 1=1 

is well-defined, closed, and bounded from below. The operator H can then be 
defined as the corresponding (unique) self-adjoint operator. It satisfies 

N N 

/\H\R^-C'') C V{H) C Q{H) = /\i7i/2(K3.c9)^ 

i=l i=l 

For Za < 1/2, V{H) = A^L^H^ (R^iC^). AU this follows from (the statements 
and proofs of) [TF', Theorem X.17] and [Et Theorem VIII.15]. See [H] for further 
references on H. We shall not have anything further to say on H in this paper, 
however, but will only study the Hartree-Fock problem mentioned above. We now 
discuss this in more detail. 

It is convenient to use the one-to-one correspondence between Slater determi- 
nants and projections onto finite dimensional subspaces of L^(R^;C*). Indeed, if 
^ is given by ^ with {u^}f^l C iJ^/^j^M^. C«), orthonormal in L'^ {R^ ; C^) , and 7 
is the projection onto the subspace spanned by ui, . . . , ujv, then the kernel of 7 is 
given by 

N 

7(x, cr; y, t) = ^ Uj (x, a)uj (y, r) . (10) 
Let p-y G L^(R'^) denote the 1-particlc density associated to 7 given by 

q q N 

(7—1 cr— 1 j — 1 

Then the energy expectation of ^ depends only on 7, more precisely, 

q(vI/,vI/) = (*,ijvl;) =5HF(^)^ 

where is the Hartree-Fock energy functional defined by 

£HF(7) = a-i{ Tr[i?(p)7] - a'' Tt[^\ - Tt[V^]] + V{^) - £x{^) . (11) 

Here, 

N N r ( ) 

Tr[i?(p)7] ^ e^-, , Tr[y7] ^ 0^-, zi,] = Za / ^ dx , 
15(7) is the direct Coulomb energy, 

^ JR3 JR3 |x - y| 

and £x{'y) is the exchange Coulomb energy, 

2^^i7r3Jr3 |x-y| 
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This way, 

(TV, Z, a) = inf { (7) I 7 e } , (13) 
V = {j: L2(M^; C«) ^ ^^(M^; C«) | 7 projection onto span{Mi, . . . , un}, 

(Notice that if one of the orbitals ui of 7 is not in H'^/'^{M?; C«), then £"^(7) = +00 
(since Za < 2/7r).) 

We now extend the definition of the Hartree-Fock energy functional in order 
to turn the minimization problem (jl3p (that is, ([3|)) into a convex problem. 

A density matrix 7 : L'^{R^\<C^) L?'{R^\'C) is a self-adjoint trace class oper- 
ator that satisfies the operator inequality < 7 < Id. A density matrix 7 has the 
integral kernel 

7(x, cr; y, t) = X! -^J^J' ^) ' 

i 

where Aj, Uj are the eigenvalues and corresponding eigenfunctions of 7. We choose 
the Mj-'s to be orthonormal in L'^{R^\C). As before, let G L^(R^) denote the 
1-particle density associated to 7 given by 

p,(x)^^^A,K(x,a)r (15) 
0=1 j 

Define 

:= {7 density matrix | Tr [-B(p)7] < +00 } , (16) 
where, by definition, for 7 written as in (|14p . 

Tr[i?(p)7] :=^A,e[u„u,]. (17) 

Notice that if 7 G then all the terms in £^^(7) (see (dJ)) are finite. Indeed, for 
7 G .4 and written as in (fTH) . 

Tr[F7] := V A,o[u„ u,] = Za / ^^1^ dx (18) 

is finite, due to ([S]). In particular, 

Uj G i/i/2(M^; C«) C L3(M^; C«) , (19) 

the last inclusion by Sobolcv's inequality [HI Theorem 8.4]. 
On the other hand, if 7 G .4 then 

G Li(M3)nL4/3(R3). (20) 

This follows from Daubechies' inequality, see [3 pp. 519-520]. By Holder's inquality, 
p-f G L^/^{^^). The Hardy-Littlewood-Sobolev inequality [H Theorem 4.3] then 
implies that I?(7) (see (fT^ ) is finite. Finally, £x{'^) < ^'(7), since 

2?(7)-fx(7) 

^i-'\ftJ^JM^ |x-y| 

Therefore, defined by (jlip extends to 7 G .4. This way, with /iq defined as 
in dS]), 

Tr[/io7] = Tr[S(p)7] - a'' Ti[j] - TT[V-f] , 
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and so 

£«F(7) = Tr[/io7] + T^il) - £x{^) , 7 e ^ . (21) 
Consider 7 G -4 and define, with as in (llSp . 

i?,(x):= / A^rfy. (22) 

|x " y| 

We have that 

i?-, e L°°(M^)nL3(M3). (23) 

This follows from © (for L°°), and (HO]) and the weak Young inequality O p. 107] 
(for L^). Next, define the operator with integral kernel 

I \ 7(x,g;y,T) 

/s:T,(x,cr;y,T) — |^r~^yj — ' ^ ' 

The operator if^ is Hilbert-Schmidt; we prove this fact in Lemma [5] below. 

Note that, using (fT4|) and the Cauchy-Schwarz inequality, {u,R^u) > {u,K~fU) 
(multiplication by R-y is in each component). Denote by bj the (non- negative) 
quadratic form given by 

b.y[u,v] -.^^ a{u,R^v) - a{u,K^v) for u, w G ^^/^(R^; C?) . 

Then, using [u.K^u) > and ([5]), 

G<bJu,u]<a{u,R^u)^aY I I ^^(y)l"(^' d^dy < a- Tr[7] e[u,u] . 

f-{ Jr3 7r3 |x - y| TT 

Therefore (by the statements and proofs of [HI Theorem X.17] and [T71 Theo- 
rem VIII. 15]), there exists a unique self-adjoint operator hj (called the Hartree- 
Fock operator associated to 7), which is bounded below (by — a^^), with quadratic 
form domain _ff^/'^(R'^; C) and such that 

{u,h^v) ^t[u,v]-o[u,v] + bj[u,v] for u,v e H^/^{R^;C'') . (25) 

The operator has infinitely many eigenvalues in [~a^^,Q) (when N < Z), and 
acss{hj) = [0, 00); both of these facts will be proved in Lemma O below. 
The main result of this paper is the following theorem. 

Theorem 1. Let Za < 2/tt, and let N > 2 be a positive integer such that N < 
Z+1. 

Then there exists an N -dimensional projection 7^^ = 'j^^{N,Z,a) minimizing 
the Hartree-Fock energy functional £ given by dm), that is, E^^{N,Z,a) m ^ 
( and therefore, in ^) is attained. In fact, 

£:HF(^HF^ ^ i?HF(^^ a) = inf {f HF(7) I 7 e ^, 7' = 7, Tr^ = N} 

= inf {£"^(7) I 7 e ^, Tr[7] ^ N} 
= inf {£"^(7) I 7 e y^, Tr[7] < N} . (26) 

Moreover, one can write 

N 

7™(x,CT;y,r) = ^(p,{x,a)(p^{y,T) , (27) 

with (fii S i?^/^(R'^; C?), I = l,...,iV, ortnonormal, such that the Hartree-Fock 
orbitals {ipi}fLi satisfy: 

(i) With h^uF as defined in (pS)) . 

h^Hwipi = Eiipi , I = 1, . . . , iV , (28) 
wit/I 0>£Ar>...>ei> — a^^ the N lowest eigenvalues of H^hf. 
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(ii) Fori^l,...,N, 

(p,; € C°°(K3\{0};C«). (29) 

(iii) For all R > and j3 < v^^^ := —eN{'2ot^^ +Sn), there exists C = 
C{R, /3) > such that for i ^ 1, . . . , N , 

|(p,(x)| < Ce^'^l'^l for |x|>i?. (30) 

Remark 1. 

(i) In fact, we prove that ()29p holds for any eigenfunction ip of h^uF , and 
pop for those corresponding to negative eigenvalues e. More precisely, if 
h^awip — eip for some e £ [sn, 0), then pop holds for (p for all (3 < := 
v/-e(2a-i +e) for some C = C{R, (3) > 0. 

(ii) Note that, in general, eigenf unctions of H^hf can be unbounded at x = 0; 
therefore p9p and pOp can only be expected to hold away from the origin. 

(iii) Both the regularity and the exponential decay above are similar to the results 
in the non-relativistic case (i.e., for the operator in see [13 ). However, 
the proof of Theorem [7] is considerably more complicated due to, on one 
hand, the non-locality of the kinetic energy operator -E(p), and, on the 
other hand, the fact that the Hartree-Fock operator /i^hf is only given as a 
form sum for Za G [1/2, 2/7r). 

(iv) We show the existence of the Hartree-Fock minimizer by solving the min- 
imization problem on the set of density matrices. This method was intro- 
duced in 23J. The same method was used in 4 in the Dirac-Fock case. 

(v) As mentioned earlier, we have to assume that Za < 2/7r; the reason is that 
our proof that Tr[i?(p)7„] is uniformly bounded for a minimizing sequence 
{7n}riGN docs not work in the critical case Za = 2/?:. 

(vi) For simplicity of notation, we give the proof of Theorem [7] only in the 
spinless case. It will be obvious that the proof also works in the general 
case. 

(vii) As will be clear from the proofs, the statements of Theorem]^ (appropriately 
modified) also hold for molecules. More explicitely, for a molecule with K 
nuclei of charges Zi, . . . , Zk , fixed at Ri, . . . , Rk G K"^, replace 'o in ^ by 

K 

0[it,w] :=^(Vfe/\,Vfc/^w) for u,wei?i/2(R3;C«), (31) 

with Vfc(x) = Zfctt/lx - Rk\,Zka < 2/n. Then, for N < 1 + X]fc=i ^fc^ 
there exists a Hartree-Fock minimizer, and the corresponding Hartree-Fock 
orbitals have the regularity and decay properties as stated in Theorem [II 
away from each nucleus. 

2. Proof of Theorem [T] 

2.1. Existence of the Hartree-Fock minimizer. The proof of the existence of 
an A'^-dimensional projection minimizing the equahties in ()26p . and that 
the corresponding Hartree-Fock orbitals {(pi}fLi solve the Hartree-Fock equations 
([28l) , will be a consequence of the following two lemmas. 

Lemma 1. Let Za < 2/tt and N eN. Then 

E^^'iN, Z, a) inf {£^^(7) | 7 e A Tr[7] < N} 

is attained. 
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Lemma 2. Let 7 S yt. Then the operator K^, defined by (|24p . is Hilbert-Schmidt. 
If Za < 2/7r then the operator h^, defined in (|25p . satisfies acsa{h~j) = [0,oo). // 
furthermore Tr[7] < Z, then h-y has infinitely many eigenvalues in [~a^^,0). 

Before proving these two lemmas, we use them to prove the parts of Theorem [T] 
mentioned above. 

Proof. For computational reasons we first state and prove a lemma in the spirit of 
[31 Lemma 1]. 

Lemma 3. Let ^ £ A, ui,U2 G H^/^{M.'^), and let 61,62 G K 6e such that 7 given 
by 

7(x,y) :=7(x,y) +7«(x,y), (32) 



7„(x,y) := 7„j^„2(x,y) = eiui(x)ui(y) + e2U2(x)u2(y) (33) 

is again an element of A. 
Then we have that 

£HFj-~-j ^ £^^(7) + a~^ei{ui,h^ui) + a^^t2{u2, hyU2) + eie2i?M , (34) 
where is given in (PS]) . and 

Ru := = U f \M^)My)-u2i^)u,iy)\^ (33^ 

2 Jb3 Jr3 |x - y| 

Proof of Lemma [H: We have that 

£:«^(7) = 5«"(7) + a-^Tr[/.o7j+ / / ^^^^^ dxdy 

JR3 Jr3 |x - y| 

7(x,y)7.(x,y)^^^^^l f f (x)p, (y) ^^^^ 



|x - y| 2 7ji,3 7r3 |x - y 



Iff 7„(x,y)7„(x,y) ^^^^ 



2 7r3 Jr3 |x - y| 
= £™{-f) + a^^ei{ui,hyUi) + a^"^ €2(1*2, ^17^2) (36) 



1 f f / / 7.(xy)7n(x,y)^^^y^ 



|x - y| 2 7r3 7r3 |x - y| 

Using dSSl), that ^^^(x) = ei|ui(x)p +e2|u2(x)P, and ([35]), we obtain ([S!]). □ 
By Lemma [T] a minimizer 7^^ e yt, with Tr[7J^^] < N, exists. We may write 

7«^(x,y)=^A,^,(x);^, (37) 

k 

with 1 > Ai > • • • > and {ipk}k C if^/^(M^) an orthonormal (in L^(]R^)) system 
(it might be finite). Extend {ipk}k to an orthonormal basis {ipk}k U {w^j^eN for 
L2(M3), with ui G ffi/2(M3). 

Let X + 1 be the first index such that Xk+i < 1. Fix j E {f , . . . , K}, choose 
u G {<Pfc}fc>/f+i U {ue}i^fi, and consider, for e to be chosen, 

7p)(x,y) := ^ Afc^fe(x)(^^(y) + ^ {^M) + ( V^jXy) + £^(y) ) ■ 

Choosing m > 1 assures that Tr[7p^] < N. Then < 7^'^ < Id for |e| small enough 
(depending on u). Since 7^^ minimizes and 7q"'-' = 7^^, 



o=|(f«^)(7.(^^: 



e=0 
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Repeating the computation for iu we get that (u, h^awipj) = 0, from which it follows 
that hjHF maps span{c/?i, . . . , ipk} into itself. Diagonalising the restriction of /i^hf 
to spanjtpi, . . . , ipk}, we can choose t/Ji, . . . , ipx to be eigenf unctions of h^aw with 
eigenvalues £„j , ■ • ■ , , n-j Cz N (numbering the eigenvalues of H^hf in increasing 
order, —a~^ < £i < £2 < • • ■)• Since Ai = • • • = = 1, this does not change ([57)1 . 

To show that, for j > K , ipj is also an eigenfunction of /i^hf (corresponding to an 
eigenvalue e„j. ) one repeats the argument above, with u £ {^k}k^i,...,K,j U {M^j^gN, 
and 

7(j')(x,y) = ^Afc</jfc(x)(^fe(y) + ^ {^j (x) + eu{x)) ( tp^ (y) + eu{y) ) . 

Moreover, the eigenvalues (of /i-^hf) corresponding to the eigenfunctions tpk 
are non-positive. In fact, if > 0, then we could lower the energy: Define 

7(x,y) = 7"^(x,y) - Xkfkix)Vk{y), then, using Lcmma[31 we get that £"^(7) = 
£:HF(7«^)-a-iAfe£„, <£HF(^HF), 

It remains to show that Tr[7^^] = N, that 7^^ is a projection, and that the 
{v^j}jLi are eigenfunctions corresponding to the lowest (negative) eigenvalues of 
hryHF (that is, to El < £2 < • • • < Eat < 0). 

Consider first the case N < Z. Assume, for contradiction, that Tr[7^^] < N. 
Let iiT S N be the multiplicity of the eigenvalue 1 in ([57|) . Since (by Lemma H]) , for 
N < Z , h^nw has infinitely many eigenvalues in [— a^^, 0) we can find a (normalized) 
eigenfunction u, corresponding to a negative eigenvalue of H^hf, and orthogonal to 
ipi, . . . ,(pK- Let e > be sufficiently small that 7(x, y) :— 7^^(x, y) + eu(x)'u(y) 
defines a density matrix satisfying Tr[7] < TV. By Lemma |3] (with ui = u, ei = e 
and 62 — 0) we get that 

£HF(^) ^ £HF(^HF) ^ ^^^^^^ < £HF(^HF) ^ (gg) 

leading to a contradiction. Hence, Tr[7^^] — N. That 7^^ is a projection follows 
from Lieb's Variational Principle (see [TT]) which we prove for completeness. If 
this is not the case, there exist indices p,q such that < Ap,Ag < 1. Consider 
7(x,y) := 7"^(x,y) + e(fq{x)ipq{y) - e(fp{x)Lpp{y) with e such that < 7 < Id. 
Choose e > if e„ < Sn and e < otherwise. By Lemma |31 we get that 

f HF(^) < £HF(^HF); 

Consider now the case Z < N < Z + 1 (and N > 2), so that N - 1 < Z. Let 
In^i denote the density matrix where 

inf I 7 e ^, Tr[7] < N ~ 1} 

is attained. By the above, Tr[7]J'^^] = — 1 and 7^^^ is a projection, so its 
integral kernel is given by 

N-l 

7"-i(x,y) = M^)<t>i{y) , 

1=1 

where the (j)i^s are eigenfunctions of H^hf ^. 
We first prove that 

inf {£"^(7) I 7 e A Tr[7] < N} (39) 

is not attained at the density matrix 7^^'Li by constructing a density matrix 7 with 
Tr[7] < N such that £^^(7) < £^^{j^^-^). Indeed, since H^hf ^ has infinitely many 
strictly negative eigenvalues (by Lemma HJ iV — 1 < Z) there exists a (normalized) 
eigenfunction u of h^HF ^ corresponding to a negative eigenvalue, and orthogonal to 
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spanj^i, . . . , (j)N-i}- Let 7 be defined by 

7(x, y) = 7"^ 1 (x, y) + w(x)u(y) . 
Then Tr[7] = N and, by a computation like in ([55]) . 

cHF/~^ cHF/ HE \ , ~1 / i \ ^ fHF/ HE \ 

Hence, 

inf \jeA, Tr[7] < iV} < inf {£"^(7) \ieA, Tt[j] < - l} . (40) 

Let 7jv be a density matrix where p9p is attained (the existence of such a 
minimizer follows, as before, from Lemma [T]). By the above it follows that — 1 < 
Tr[7Ar] < N. We now show that there exists a minimizer 7^^ with Tr[7^^] = N. 

The integral kernel of is given by 

7A,(x,y) = ^ Aj(pj(x)(^j(y) , 
j 

where 1 > Ai > • • • > and the (pj^s are (orthonormal) eigenfunctions of /i-y„. If 
Tr[7Ar] < N we can define a new density matrix 7 with Tr[7] < N and (7) < 
£^^{jiy). Indeed, if Tr[7jv] < N (and bigger than — 1) then there exists a (first) 
jo such that < Aj„ < 1. We define 7 with integral kernel 

7(x,y) =7Ar(x,y) +ripj„{x)tpj„{y) , (41) 

with r = min{l — Xjg,N — Tr[7jv]} > 0. Recall that h^j^ipj = En^ipj, < 0, for 
all j. By Lemma [3] we have that 

£«^(7)=f«^(7A^)+a-V.„,„. 

If Euj^ < 0, it follows that £^^(7) < £^^(7Ar). On the other hand, if £„^^ = 0, then 
fHF(^) = £.HF(^^-)^ g^^^ Tr[7Ar] < Tr[7] < N. Either Tr[7] = N, in which case we 
let 7^^ := 7, and, as above, we are done. Or, we repeat all of the above argument 
on 

jo 

7(x,y) = W'^i(y) + ^j'^iW'^i(y) • 

J = l j>3o 

Since the trace stays bounded by N, this procedure has to stop eventually. Hence, 
with 7^^ the resulting density matrix, Tr[7^^] = N and by Lieb's Variational 
Principle it follows (as above) that 7^^ is a projection. 

Finally, let {(pj} be the eigenfunctions of ^..^hf, now numbered corresponding to 
the eigenvalues £1 < £2 < ■ • • , where £1 is the lowest eigenvalue of h^HF. We know 
that, for some ji, . . . , jat e N, 

N 
k=l 

Suppose for contradiction that {£ji, . . . , £j„} 7^ {^ii • ■ • i^w}- Then there exists a 
fc e {1, . . . , iV} with £jj^ > £fc. For 5 e (0, 1) define 

7(x,y) = 7"''(x,y) +(5(^fc(x)v3fe(y) - (5<^j Jx)(^j, (y) . 

By Lemma [31 

£HF(~) ^ ^HE(^HE) ^ _ _ ^2^^^^^^^ < ^HF(^HF) ^ 

where the last inequality follows by choosing 5 small enough. 
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It remains to prove that £i, . . . , Eat are strictly negative. For N < Z this foUows 
directly from Lemma [H In the case Z < N < Z + 1, assume, for contradiction, 
that Eat = 0; then the density matrix 

7(x,y) 7"^(x,y) - (^Ar(x)(^Ar(y) 

satisfies £^^(7) = £^^(7^^) (by Lemma [3]) and Tr[7] iV — 1. This is a contra- 
diction to (001). 

This finishes the proof of the first part of Theorem [TJ □ 
It remains to prove Lemma [T] and Lemma [21 

Proof of Lemma]^: We minimize on density matrices following the method in [23]. 
In the pseudorelativistic context one faces the problem that the Coulomb potential 
is not relatively compact with respect to the kinetic energy. This problem has been 
adressed in [3] and we follow the idea therein. 

The quantity E^^{N, Z, a) is finite since for any density matrix 7, with Tr[7] < 

£™(7) > a"i{Tr[£;(p)7] ~ a-~^ N ~ Tt:[V -f]] > -a-'^N. 

Here we used that 15(7) - £x{-f) > 0, and (0) (see also ^ and JTH])). 

Let {7n}J^Li be a minimizing sequence for E^^{N, Z, a), more precisely, 7„ G ^ 
(with A as defined in (HH)), Tr[7„] < N, and 5"^(7„) < E^^{N, Z, a) + 1/n. 

The sequence Tr[i?(p)7„] is uniformly bounded. Indeed, for every n e N, using 
®, 

E^''iN,Z,a) + 1 > £"^(7n) > Tr[£;(p)7„] - a'^N - TT[Vjn]} 
> a-\l - Za^) Tr[i;(p)7„] - a'^N . 

The claim follows since Za < 2/tt. It is this argument that prevents us from proving 
Theorem [T] for the critical case Za = 2/7r. 

Define 7„ := £:(p)i/2^„£:(p)i/2. Then, by the above, is a sequence of 

Hilbert-Schmidt operators with uniformly bounded Hilbert-Schmidt norm. Hence, 
by Banach-Alaoglu's theorem, there exist a subsequence, which we denote again 
by 7„, and a Hilbert-Schmidt operator 7(00)1 such that for every Hilbert-Schmidt 
operator W, 

Tr[T^7„] ^ Tr[iy7(o,)] , 71 ^ 00 . 

Let 7(oo) •= -E'(p)~^^^7(oo)-K(p)""'^^^- We are going to show that 7(oo) is a minimizer 
of (in fact, of a£^^ , which is equivalent). We first prove that 7(oo) G then 
that is weak lower semicontinuous on A. 

Let {TZ-felfeeN be a basis of L'^{R^) with ipk e H^^'^iM.^). Then, for all fc e N, 

lim (V'fe,7n^fe) = lim {^Pk, E{p)-^^^nE{p)-^^^^Pk) 

n—*oo n—'oo 

From this follows, by Fatou's lemma, that 

Tr[7(oo)] = y](V'fe,7(oo)'/'fe) < liminf V(?/'fc,7„V'fe) = liminf Tr[7„] < N, 

k k 

and 

TT[E{p)^/\^-,E{py/'] < liminfTr[£;(p)i/27„i?(p)i/2] < ^. 

n — ►oo 

Since also < 7(00) < Id we see that 7(00) € A. 

To reach the claim it remains to show the weak lower semicontinuity of the 
functional . As mentioned in the introduction, the spectrum of the one-particle 
operator ho, defined in Q, is discrete in [— a^^, 0) and purely absolutely continuous 
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in [0,00). Let A_(a) denote the projection on the pure point spectrum of /iq and 
A+(a) := Id - A_(a). We write 

a£"P(7„) = ri(7„) + T2{in) + aTsi^n) , (42) 

with 

Ti(7n) = Tr[A+(a)/ioA+(a)7n] , T2{-fn) - Tr[A_(a)/ioA_(a)7«] , 

We consider these three terms separately. 

For the first term in ([HI), fix (as above) a basis {V'fcjfeeN of L^(R'^), with 
{V-felfceN C ifi/2(R3), Defining 

Jk := (A+(a)/ioA+(a))^^Vfc, 

we have that 

Ti(7„) = Tr [(A+(a)/^oA+(a))'/'7„(A+(a)/ioA+(a))'/'] 
= E(/fe'^"/'=) = E(^(P)"'^Vfc,7n£;(p)-^/Vfe)- 

k k 

Since the projection 

is a non-negative Hilbcrt-Schmidt operator, we find, by Fatou's lemma, that 
liminf Ti(7„) = liminf Tr[iJfe7„] > Tr[i7fc7(oo)] = ^'1(7(00)) ■ 



77 — >00 

A; fe 



As for the second term in (|42p . we have linx^^oo r2(77i) = 12(7(00)) since the 
operator K-{a)hQK^{a) is Hilbert-Schmidt; see Lemma[7]in Appendix Kl 

Finally, for the last term in ((1^ . following the reasoning in ^ pp. 142-143] (here 
we need that TV e N), we get that 

liminf T3(7,7) > T■i{-f,r^^) . 

n — ^00 

This finishes the proof of Lemma [U □ 



Proof of Lemma[E- In order to prove that is Hilbert-Schmidt it is enough to 
prove that its integral kernel belongs to L^(R^). We have that (see (fM|) and (fH|) ) 

/ |i^^(x,y)|2dxdy= / M^ElZlf rfxdy (43) 
Jro Jr6 |x - y| 

^ ^ r ^. (x)..(y)^ ^ _ 

The last integral can be estimated using the Hardy-Littlewood-Sobolev, Holder, 
and Sobolev inequalities (in that order), to get 

Ij,k < \\ukUj\\y2 < hfcllall^jlls < C'll"fc|lHi/2||wj||^i/2 . (44) 
Inserting (jU]) in we obtain (since E A) 

[ \K^i^,y)f dxdy < C ^ •^j^fell^*fellHi/2 hi 11^1/2 = C XjWu^f^./^) 

= C(Tr[i?(p)7])' <oo. 

To prove the statement on the essential spectrum, define + aK^. Since 

is Hilbert-Schmidt, and CTcss(^o) — [0,oo) (see the introduction), it is enough 
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to prove that {h^ +77) ^ — {ho + rj) ^ is compact for some r] > large enough ^Hl 
Theorem XIII.14]. Since V{hQ) = V{h^) C V{R-y), we have that 

{hj + ri)-^ - {ho + 77)^1 = -~{h-i + rjy^aR^{ho + r;)"^ . (45) 

From Tiktopoulos's formula (see [22l (II. 8), Section II.3]), it follows that 

{ho + vr' 

= (r(p) + ry)-i/2[i _ (T(p) + fjr'/^V{T(p) + ry)-!/^]-! (r(p) + r^y'/^ . (46) 

Since, by jS]), ||(T(p) + r/)-i/2^i/2|| ^ f^j. ^-q, < 2/n and 77 > the right side 
of ((46|) is well defined. Inserting (|46| in (|45|) one sees that it suffices to prove that 
Rj{T{p) + r])~^/^ is compact. That this is indeed the case follows by using [T^ 
Theorem XI. 20] together with the observation that, for e > and 77 > a"^, Rj and 
(T(p) + 77)^1/2 (as a function of p) belong to the space L'^+'^(K''^) (for R-^.. see (P5| ). 

Finally, we show that if Tr[7] = < Z then hj has infinitely many eigenvalues 
in [— Q;~"'^,0). By the min-max principle 20. Theorem XIII. 1] and since acssih^-f) = 
[0,00), it is sufficient to show that for every n G N wc can find n orthogonal 
functions ui, . . . , u„ in L'^iM.^) such that {ui, h^Ui) < for i = 1, . . . ,n. 

Let n E N. Fix S 1 — N/Z and let ho^s be the unique self-adjoint operator 
whose quadratic form domain is H^/^CR"^) such that 

{u, hQ^sv) = t[u, v]~6 v[u, v] for u,v e i7^/^(R^) . 

By [71 Theorems 2.2 and 2.3], (Jcssiho^s) = [0, 00). Moreover, ho s has infinitely 
many eigenvalues in [— q;~^,0). This follows by the min-max principle and the 
inequality ho^s < a/2{—A) — 5Za/\x\. Hence, we can find ui, . . . ,u„ spherically 
symmetric and orthonormal such that (u^, ho^sUi) < for i = I, . . . , n. Then, by 
the positivity of K~^, by Newton's Theorem [12l p. 249], and since Tr[7] = N we 
get, for i — 1, . . . ,n, that 

{ui,h^Ui) < t[ui,Ui] - o[ui,Uj] -I- a{ui,RjUi) 

N 

< t[ui,Ui] - v[ui,Ui\ + — t)[ui, Ui] = (Ui, ho.sui) < . 
The claim follows. □ 



2.2. Regularity of the Hartree-Fock orbitals. Here we prove that any eigen- 
function of /i^hf is in C°°(MP \ {0}). 

Proof. Let cp he a solution of h^HFip — eip for some e e K. Then (p belongs to the 
domain of the operator and in particular to iJ^/^(M^; C). We are going to prove 
that ip e H''{n) for all bounded smooth C \ {0} and aU fc g N. The claim 
will then follow from the Sobolev imbedding theorem [2, Theorem 4.12]. We will 
use results on pseudodifferential operators; see Appendix |B] We briefly summarize 
these here. 

1) For all fc,£ e R, E{pY maps H''{R^) to H''-'^(R^). 

2) For all k,£ € R, and any x G Cg°(R^), the commutator [x,E{pY] maps 
H'^iR^) to H'^-'^+^R^). 

3) For all k,£,m e R and Xi>X2 & C^(R) with suppxi n suppx2 = 0, 
XiE{pYx2 maps H''{R^) to H"^{M.^). Such an operator is called 'smooth- 
ing'. 

Fix a bounded smooth subset of R"^ \ {0}. We proceed by induction on fc S N. 
Assume that ip e H''{n) for some fc > 0, i.e., XV' e -^''(R^) for aU x e C^{^)- 
Notice that H''{R^) = D{E{pf). 
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Since G H''+^{M?) is equivalent to x'P G D{E{-pf+'^), and D{E{-p)''+'^) = 
D[[E{tpY+^)*), it is sufficient to prove that G D{{E{Tpi)^+^)*), or cquivalently, 
that there exists v £ L^(R'^) such that 

(X^, i?(p)'^+V) = (^', /) for aU / e H'=+1(M3) . 

Let / e i7'=+i(R3). Then 

(X^,£;(p)'^+V) = e(^,i?(p)-ixi?(p)'+V) 

= (£ + E{p)'\Eip)'^+'f) + E{p)-\Eipf+'f ) 

-b.,H.{^,E{p)-'xEip)'+'f), (47) 

where we use that h^uvip — e^p. We study the terms in (|47l) separately. In the 
following, X denotes a function in C^(il) with x = 1 on suppx- 
For the first term in (|T7l) we find that 

{if,E{p)-\E{pf+^f) = (xi?(p)-V,i?(p)'+\/) 

= ([x, E{p)-']ip, E{p)'^+'f) + iE{p)-'xV>, E{pf+'f) . (48) 

Since xv G H^{M.^) by the induction hypothesis, we have that E{p)^^x^ G 
^fc+ij-jjs-) g^j^j hence there exists wi G L^(M^) such that 

{E{p)-\v,E{pf+'f)^{w^J). 
It remains to study the first term in (|48l) . We have that 

([x,i?(p)-^]^,i?(p)'^+V) 

= ([x, i?(p)-']x^, i?(p)'+V) + ([x, i?(p)-'](i - x)^, i?(p)'+V) . 

Since xf G i/'^(]R"^) by the induction hypothesis, it follows from Proposition [2] that 
[x, £^(p)~^]X¥' belongs to i/''+^(M'^). On the other hand since the supports of x 
and X are disjoint the operator [x, i?(p)^^](l — x) is a. smoothing operator. Hence 
there exists a. ■W2 £ Lp'iM?) such that 

{[x.E{p)-']^,E{pf+'f)^{w2j). 

As for the second term in (|T7)) . we find, with x as before, 

X>{^,E{p)-\E{pf+'f) ^ {^,VE{pr^xE{pt+'f) 

= {x^,VE{p)-\E{p)''+'f) (49) 

+ ((1-X)^, VE{p)-\E(pf+^f). 

Since x bas support away from zero, Vxip G H^'iM?) and hence there exists G 
L2(M3) such that 

(X^, VE{p)-\E{pf+'f) = {W3j) . 
For the second term in (|49p we proceed via an approximation. Let {</?n}^i C 
Cg°(R^) such that ip„ ^ ip,n ^ oo, in ^2(^3)^ gj^^^g _ x)^^^(p)"^x£^(p)''+ V 
belongs to L'^{M.'^), we have that 

{^,{l-x}VE{py\E{p)''+'f)= hm ((^„,(l-x)V^i?(p)-'xi^(p)'+V). 

n—^-\-oo 

For each n G N, V{1 - x)<fn & iJ™(M^) for all m, since G C^(M^), and V maps 
£f'=(R3) into ir'=-i(R3) for all k. Therefore, E{p)''+^xEip)~'^V{l-x)^n S ^^(R^), 
and so 

(<^„,(i-x)^^i?(p)-'xi?(p)'+V) 
= (£;(p)^+ixi?(p)-V(i-x)^„,/) 

= iEip)'^+'xEip)-\l - x)Eip)E{p)-'V^^, f) . 
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Here E{p)~^V is bounded by and x£'(p)~^(l — x) is a smoothing operator by 
the choice of the supports of x and %. It then follows that {£'(p)'^+^xi?(p)^^(l — 
X)E{p)E{p)~^Vipn}n<£N is a uniformly bounded sequence in L'^{R^) and hence 
there exists W4 e L^(R'^) such that 

hm ((p„, ((1 - x)VE{p)-\E{pf+'f) = K, /) . 

n — *-\-oo 

For the third term in (IT71) . we have to separate the cases k = and fc > 1. 
Let fc = 0. The terms R^hflp and K^HFip belong to L^(M^), since R^hf e L°°(M^) 
(see and X^hf is Hilbert-Schmidt (see Lemma and therefore 

b^^F{^,Eip)-'xEiP)f) = a iEip)xEip)-\R^^F - K^^f)^, f) . 

Assume now fc > 1. With x as before, 

b^HF{^,E{p)-\E{p)''+'f) = a{x{R^HF-K^HF)^,E{p)-\E{p)''+'f) (50) 

+ a ((1 - x){Rj^^ - K^^^)V, Eip)-\E{p)''+\f) . 

By the induction hypothesis and Lemma [6] (see Appendix E| we have that xR-y^^V 
and xK-yFiF ip belong to H''{M.^). Therefore there exists W5 e L'^{M.^) such that 

(x(i?^HP ~ K^hf)^, Eip)-'xEip)''+'f) = (W5, /) . 
For the second term in (jSOp we find, since R^nFcp, K^kflp S L^(M'^), that 

((1 - x)(i?7H- - i^7HP)<^, i?(p)-ixi^(p)'+V) 

= (Xi?(p)-l(l - X)(i?7HP - X^HP)^, i?(p)'= + V) , 

and the result follows since x-E'(p)^^(l — x) is a smoothing operator. □ 

2.3. Exponential decay of the Hartree-Fock orbitals. The pointwise expo- 
nential decay ([30|l will be a consequence of Proposition [T] and Lemma |4] below. 

Proposition 1. Let 7^^ he a Hartree-Fock minimizer, let H^hf be the correspond- 
ing Hartree-Fock operator as defined in (|25p . and let {fiJfLi be the Hartree-Fock 
orbitals, such that 

h^HFifi = e.iipi , I = 1, . . . , TV , 
with > En > • • • > El > — ot^^ the N lowest eigenvalues of h^uF . 



(i) Let :— eAr(2Q; ^ + En)- Then ipi G P(e^l ' ') for every (3 < v^n '^"■^ 
ie{l,...,N}. 

(ii) Assume h^HFf — eip for some e £ [£Ar,0), and let '■— \J —si^a^^ + e). 
Then (p G I?(e^' ' ') for every P < v^. 



Lemma 4. Let E < and ve := ^1 - E{2a~^ + E)\ = yJ\or'^ - {E + ar'^Y\ 
Then the operator T(— iV) — 
tegral kernel of its inverse is l 

(T-i?)-i(x,y) =G£(x-y 



Then the operator T(— iV) — E = \/ —/S. + a — a. ^ — E is invertible and the 
integral kernel of its inverse is given by 

{E + a-'^)e-'"^\''-y\ Xi(a-i|x-y|) 



a 



-1 



47r|x-y| 27r2 |x-y| 



where Ki is a modified Bessel function of the second kind [I]. 
Moreover, 

n<r (.\<r I «'^ ^i(^'^|x|) 

< G^;(x) < C.,^,^^ + ^ , (52) 

e'^l-lGis G L«(R^) /ora?/ /3 < j^b and gG [1,3/2). (53) 
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Proof of Lemma^: The formula ((5T|) for the kernel of {T — E) ^ can be found in 
m eq. (35)]. 

The estimate (|52|) is a consequence of the bound 

j^i(a^i|-|) ^ e-^N e--^l-l 

— n — ^ - ^"'"^^^ • 

This estimate, on the other hand, follows from Newton's theorem (see e. g. |12)). 



|x-y| 47r|y| 



■dy 



Jr3 |x-y| 47r|y| 47r|x| J^s |z| 

The last integral is finite since ve < oT^ , using the following properties of K\ (see 
[g 8.446, 8.451.6]): 

Kx{t) < -i- for aU t > 0, (54) 
r I 

and for every r > there exists c,- such that 

i^ift) < c,.^ for all t > r . (55) 
The estimate ((55)) is a consequence of (|5^ . and (|55)) . □ 

Before proving Proposition [TJ we apply it, and Lemma |H to prove the pointwise 
exponential decay, i.e., the estimate in (|30p . 

Proof of Theorem [I] (iii) : Fix i £ {!,..., N}. If Za < 1/2 we can rewrite the 
Hartree-Fock equation P5|) as 

Za 



(V-A + a 2 - Q. i)(^j Eji^j + — - Q;_R^HF(pi + Q;i4:^HF(pi . (56) 

The idea of the proof is to study the elliptic regularity of the corresponding parametrix. 
By Lemma m we find that 

f Z (x 

(^i(x) = / {T ~ £Ny^{^.,y)[{ei- £N)fi + -r—r^i- aR^iiwLpi + aK^iiFLpi]{y)dy . 

In the case 1/2 < Za < 2/n, on the other hand, the operator of which we 
are studying the eigenfunctions cannot be written as a sum of operators acting on 
L'^{M.'^) and hence we cannot write directly the equation as in However, 
since the eigenfunctions are smooth away from the origin we are able to write a 
pointwise equation for a localized version of (pi. In fact, let x G C°°(R^) be such 
that < X < 1 and 

1 if |x| > 1, 
if|x|<l/2, 

and let, for R > 0, x_r(x) = x(x/i?). We will derive an equation (similar to ([K^ ) 
for T{—iV)(xRfi)- Indeed, for every u G we have that 

{u,h^iiF{xR(pz)) = c{u,XR^(Pi) - a^^{u,XRfi) - v{u,XR(Pz) + b^iiF{u,XR'Pi) 
= iXRU, h^uFcpi) + e(u, XRfi) - ^{xru, (fii) 
+ b^iiF{u,XR(Pt) - b^iiF{xRU,ipi) . 

Note that 

e(w, Xfl<P«) - i{XRU, (Pi) = {u, [E{p),XR](Pt) , 



X(x) = 
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where [E{p),xr] is a bounded operator in L^(]R'^) (see Appendix IB|) . and 

with K. the bounded operator on L^(]R^) given by the kernel 

/C(x, y) a 2^ LP J (x)(p^- (y) — — — . (57) 

Therefore there exists w E L'^{M.'^) such that 

Ku,XR'fii) = {ei + a^^){u,XR'Pi) + v{u,XR<fi) - b^iiF{u,XR'Pi) 
+ {u, [E{p),XR]ipi) + (u,IC(pi) = {u,w) . 

Hence XRfi S H^{R'^) and we can write the pointwise equation 

(V-A + q;-2 _ a'~'^)xR(Pi = eiXRfi + j^XRVi - aR n^XR-^i 

|x| 

+ aK^iiF{xRPi) + [E{p),XR]Vi + 'ICLpi . (58) 

This is the substitute for (|56p in the case 1/2 < Za < 2/7r; if Za < 1/2, the proof 
below simplifies somewhat, using (j56p directly. 
By Lemma m (155)) implies that 



Xi?(x)(^,(x) 



/ (T-en) ^{x,y)[^XR'fi - OiR-y^''XR'f i + oiK^HF{xR(pi) 

+ (si - eN)XRV>i + [Eip),XR]fi + !<^V>i\ (y) dy . (59) 
We will first show that, for all i? > and /? < i/^^, 

Xi?,<^ie^l-I e LP(M3) + L~(]R^) forpe[2,6), (60) 

and then, by a bootstrap argument, that XRVi^^^ ' G i°°(M'^), which is the claim 
of Theorem [T] (iii). 

We muhiply jSS]) by Xj?/2(x)e''l^l. Using that |(Za/|y|)xi?,(y)| < {Za)/R for aU 
y e R^ (ESI), (El, and ^ (recall (ET]), that e i/i/2(M3)^ and (©) we get, for 
some constant C — Cr^o > 0, that 

N 

|xi?(x)^.(x)e''l-l| <CxK/2(x)e'5|-l / (T - sat)-! (x, y) [|^,(y)| + ^ (y)|] dy 

+ Xi^/2(x)e''l"l / (T-£A,)-i(x,y)([£;(p),Xfl]</50(y)^y|- (61) 

We will show that the first term on the right side of (pT|) belongs to LP(M^) for 
p e [2,6), and that the second belongs to L°°(R3). This will prove (|Sn| . 
The first term on the right side of (|6ip is a sum of terms of the form 

hf{x) :-Xi?/2(x)e'5|"l j^^(T-£A.)-i(x,y)|/(y)|dy, (62) 

with / such that, by Proposition [Tl /e*^' ' ' G iy^(M^). By Lemma|l]we have, using 
g|x|-|y| < gl^-yl, that 

\hf{x)\<C f e^l-'''lG,„(x-y)e^l''l|/(y)|dy. 



From Young's inequality it follows that hf G LP{M.^) for all p g [2, 6), since /3 < i/^jv , 
so (by Proposition [I]) /e'^l ' I G i^^j^S) ^^^^ (^y LemmaS]) e^^l ' 'Ge„ £ i«(]R3) for aU 
9 G [1,3/2). 
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We now prove that the second term on the right side of ((M]) is in L°°(R'^). This 
follows from Young's inequality once we have proved that 

e'^l ■ I [E{p),XR]'P^ e LP{R^) for p G [2, oo) , (63) 



since 



e/3|x| 



(T - e^)-i(x, y){[Eip),XR]ip^) (y) dy 
< [ e^\--y\C{x-y)e^\y\\[E{p),XB]^^\iy)dy, 

and e/^l-lGe^ G L'i{R^) for q e [1,3/2). 

To prove ([50)1 it therefore remains to prove To do so, we consider a new 

localization function. Let 77 G C|5"(R'^) be such that < 77 < 1 and 



77(x) 



1 if i?/4 < |x| < 3i?/2 

if |x| < i?/8 or |x| > 2i?, 



and consider the following splitting 

e^l ■\[E{p),XR]v^ = ef'\ ■ \v[E{p), Xi?](w.) + e'" ' 'r/[i?(p), Xfl]((l - 

+ e'^l ■ I (1 _ ^) ^E{p), xr] + e^l ■ I (1 - 77) [i?(p) , Xfl] (1 - ?7)^. • (64) 

Since -qi^i G iJ*'(M'^) for all fc G N (as proved earlier), [E{p),XR\{v'-Pi) belongs to 
H^{M.^) for all A; G N. Hence, since r] has compact support away from x = 0, 
the first term on the right side of l|64p is in L^iR?) for p G [l,oo] by Sobolev's 
imbedding theorem (the term is smooth). 

For the second term in we proceed by duality: We will prove that 



^(x) (e^l-l77[i?(p),Xfl]((l-r,)^,))(x) 

defines a bounded linear functional on L^{R^) for any q G (1,2]. It then follows 
that V G LP{R^) for all p G [2, 00). 
Note that TS, 7.12 Theorem (iv)] 



47r2 



/ / g2(a-^|x-y|)dxdy for .g G ^(R^) , (65) 

Jr3 Jr3 |x - yr 



where iir2 is a modified Bessel function of the second kind (in fact, K2{t) = 
^ti[t-^Ki{t)]), satisfying [I] 

K2{t) < Cr^e-^ for < > 1 . (66) 

Let / G C^f (R"^). Using and polarization, we have that 

TWVXx) dx = (/, e^l ■ lr;[i?(p), Xi?J((l - 

3 

ff Xfl(x) - Xi?.(y) ^ -1, |x 

^lk_^^R,. |x-yP 

X [/(x)e'3|-lr7(x)(l - ry(y))^,(y) - /(y)e'3|^l77(y)(l - ry(x))v.,(x)] dxdy , 
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by the properties of x and ?7- Hence, 



/(x)V(x) dx 

( 

<Cr(( |/(x)|e'^l--''li^2(a-i|x-y|)e'^lyl|^,(y)|dxdy, 

J J|x-y|>fl,/4 

<CrII |/(x)|e'^l--''lif2(«->-y|kfl/4(|x-y|)e'5|''l|^.(y)Mxdy. (67) 



Note that, since (3 < Vf,^ < a ^, ([M]) imphes that e^^^'^K2{a ^\ ■ \)xR/i is in 
U{m!^) for all r > 1. Since (by Proposition [l]) ' ^ ip^ L'^ [M?) , Young's inequality 
therefore gives that 

{e^\ ■ \K2{a-^\ ■ \)xb.,a) * {e^^ ' ' l^^l) e L%R^) for aU s e [2, ex.) . 

This, (dH), and Holder's inequality (with l/q+l/s = l) imply that, for all / G 
Cg"(M3) and all q G (1,2] 

7MV^(x) dx < Cfi||(e^l ■ li^2(a-^| • \)XR/4) * (e^l 11/11, ■ 

By density of C^(K'^) in L*(M^), it follows that ip defines a bounded linear func- 
tional on L'?(M3) for any q € (1, 2], and therefore, that tJj £ ^^(R^) for aUp e [2, oo). 

Proceeding similarly one shows that the two remaining terms in (j64p are also in 
Lp(R3) for aUpe [2,cx)). 

This finishes the proof of and therefore of ((Sn|) . 

Finally we prove that XRfi^^^'^ ^ L°°{M.'^). We start again from (pT|) . We 
already know that the second term is in L°°(R'^). The first term is a sum of terms 
of the form (see also (p^) 

/i/(x)=Xfl/2(x)e'5W / (T-eA.)-i(x,y)|/(y)|dy, 



with / e L2(M3) and XJ?,/4e'^'''/ e ^^(R^) + L°°{R^) for p e [2,6) by what just 
proved, replacing i? by i?/4 in ([50]) . We find that 



/i/(x) < Xfly2(x) / e'3|-yl(T-£A.)-i(x,y)e''lylxK/4(y)l/(y)|dy 

+ Xi^/2(x) / e'3|-yl(T-e^)-i(x,y)e/'lyl(l-XK/4)(y)l/(y)|dy: 



and, again by Young's inequality, we see that both terms arc in L°°(R''^). Notice 
that in the second integrand |x — y| > i?/4. 

This finishes the proof of Theorem [T] (iii). □ 

It therefore remains to prove Proposition [TJ 

Proof of Proposition]^: We start by proving (i). It will be convenient to write the 
Hartrce-Fock equations h^aFipi = ei(pi, i = 1, . . . , N , (see (pS])) as a system. 
Let t be the quadratic form with domain [H^^^CR)]^ defined by 

N 

t(u,v) = ^t(u„v,) for all u,v e [H^/^{R^)f , 
1=1 

where u^ denotes the z-th component of u e [H^^^{W^)]^ and t is the quadratic 
form defined in ([7]). Similarly we define the quadratic forms v, r-^ and k-^, all with 
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domain [H^^^{R^)]^, by 

N N 
v(u,v) = ^D(ui,Vi) , r^(u,v) = Q;^(ui,i?^Vj) , ^(u, v) = q;(u, K-^v) , 

i=l 1=1 

with defined in Rj defined in ((22|) . and the TV x A^-matrix given by 

JR3 |x - y| 

The effect of writing the Hartree-Fock equations as a system is that is a (non- 
diagonal) multiplication operator. This idea was already used in [TS]. Note that 
(K^)ij G L^{R^) n L°^{R^); the argument is the same as for (f^ . 

Let finally E be the N x N matrix defined by (E)jj — ~ei6ij. 

We then define the quadratic form q by 

q(u, v) = t(u, v) - v(u, v) + r^(u, v) - k-y(u, v) + (u, Ev) . (68) 

One sees that the quadratic form domain of q is [H^^^{W^)]^ , that q is closed (since 
t is closed), and that there exists a unique selfadjoint operator H with I?(H) C 

(u, Hv) = q(u, v) for all u G [H^^'^(E?)]^ , v € X>(H) . 

Notice that the vector $ = (i^i, . . . , (pjy) satisfies H$ ~ 0. 

Let W{k), k G C^, denote the multiplication operator from a subset of [L^(]R^)]^ 
to [L-^(R'^)]^ given by /(x) i-^ e"' ''/(x). Instead of proving directly the claim of 
the proposition, we are going to prove the following statement, which implies the 
proposition: 

$ e V{W{k)) for ||Im(K)||R3 < ly,^ , (69) 

where $ = {(pi, . . . , <fN)- Here, k = Rc(k) + ilm(K) with Re(K), Im(K) € . 

We know that W{k)<^ is well defined on [L^(R^)]^ for k g and we need to 
show that it has a continuation into the 'strip' '^v^^ ■, where 

:={KeC3|||Im(K)||R3 <t}. 

We shall also need Sq-i; note that S^-i D ^i^ejv- '^^^ ^^"^^ O'Connor's 
Lemma (see Lemma [5] below). 

Starting from the quadratic form q defined in (j68p we define the following family 
of quadratic forms on [H^/'^{M?)]^ : 

q(«;)(u,u) ■.= ^{W{-K)n,W(~K)n), 

depending on the real parameter k e R'^. From the definition, 

q(K)(u, u) = t(K)(u, u) — v(u, u) + r-y(u, u) — k-y(u, u) + (u, Eu) , 

where 




One sees that qi^n) extends to a family of sectorial forms with angle < and 
that q(K) is holomorphic in the strip Sq-i (indeed, ||Im(K)||R3 < is needed to 
assure that the complex number under the square root in (|70p has non-negative 
real part for all p g R'^). Moreover, q(K) is closed. Indeed, it is sufficient to prove 
that the real part of q(K) is closed, which will follow from 

v(u,u) -f r^(u, u) +k^(u, u) -|- (u, Eu) < feRe(t(K))(u, u) + K{vL,n) , (71) 

with 6 < 1, if > and Re(t(K)) closed. We now prove ([7T|) . We already know that 

r.^(u,u) +k^(u,u) + (u,Eu) < K' {\x,Vi) for K' > Q . (72) 
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By © we find 

N „ 

v(u,u)<(Za)jE / |p||u,(p)pdp 

^ f r 

<{Za)^RY,[ \Mp)\'dp+ / |p||u,(p)pdp]. 

^ »=i J\p\<R J\p\>R 

Let (5 > be sucli tliat Zaf (1 - Sy^ < 1. 

Af » 3 

Re(t(K))(u,u) = V / |a-2+ Vfe -K,)2|'/'cos(^(p,K))|u,(p)|2dp 



(73) 



a (u, u) 



witli 



2cos^(e'(p,K)) - 1 = 



a 



ELfe-Re(«,))'-(Im(«:,))2) 



^3 



there exists R > such that cos(6'(p, k)) > (1 — (5) for |p| > R. Hence we find that 

N „ 3 
Re(t(«;))(u,u)>(l-<5)E / \a-^ +J2iP, - ^,?\'^' \Mp)\' dp 



=1 J\p\>R 



a ^(u, u) 



A' 

>{1-5)Y, (|p|-C)|u,(p)pdp-«-i(u,u), (74) 

,^i^|p|>i? 

with C > \\Ke{K)\\]i3 . The estimate in (|7ip follows combining ((7^ with (|73p and 
(El. 

The fact that Re(t(K)) is closed follows from 

/(|p|-C)|u,(p)pdp<Re(t(^))(u,u)<E /(|p|+C)|u.(p)pdp, 

with C > 2a-^ +Re(K). 

Hence, q(K) is an analytic family of forms of type (a) ([S] p. 395]). The associated 
family H(k) of sectorial operators is a holomorphic family of operators of type (B) 
and has domain in a subset of [iJ^/^(]R^)]^. 

We are interested now in locating the essential spectrum of H(k). Since 
is a Hilbert-Schmidt operator, the essential spectrum of H(k) coincides with the 
essential spectrum of the operator associated to 

t(K)(u, u) — v(u, u) + a r-y(u, u) + (u, Eu) . 

Notice that the operator associated to this quadratic form is diagonal. Proceed- 
ing as in the proof of (Tcss(^7) = [0,oo) (Lemma [5]), one sees that (Toss(H(k)) C 

cress(T(K) — En) with T(n) :— \^ ar"^ + Ej=i(Pj ^ ^lY' ^ oT^ ■ Hence we find that 



CTcss(H(k)) C {z e C I Re(z) > _ ||im(K)||23 - a ^ - Eat} . 

Hence 0, eigenvalue of H(0), remains disjoint from the essential spectrum of H(k) 
for all K € (recall that Y^y C S„-i) . 

Since H(k) is an analytic family of type (B) [501 P-20] in Ei,^, is an eigenvalue 
of H(0) and moreover, remains disjoint from the essential spectrum of H(k), it 
follows that is an eigenvalue in the pure point spectrum of H(k) for all k, G E^^^ 
(reasoning as in [501 page 187]). Let P(k) be the projection onto the eigenspace 
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corresponding to the eigenvalue of the operator H(k). Then P(k) is an analytic 
function in and for k G S,/ and Kn € R we have 

P(k + Ho) = Wino)Pin)Wi-Ko) . 

Here we used that W{—ko) is a unitary operator. The result of the lemma follows 
by applying Lemma [H] below to W{9) := e''^'^'^ with k e R'^, ||k||r3 = i/^^, and 
6* G {z € C| |Im(z)| < 1}. Notice that W{9) = VK(6'k) and that the projection 
P{0) := P{9k) is analytic and satisfies P{0 + Oq) = W{0o)P{e)W{-0o) for Oq £ K. 
This finishes the proof of (i). 

To prove (ii), we can work directly with the Hartree-Fock equation, since, from 
(i), the function K^nwip is exponentially decaying. Therefore, let 

q[u,v] = {u,h^HFv) - s{u,v) for u, u G i?^/^(R^) , (75) 

and note that, by assumption, is an eigenvalue for the corresponding operator (ip 
is an eigenfunction) . Define, for k G M'^, 

q(K)[u, w] = c\[W {—k)u,W {~k)v] 

— t{K)[u, v] — v[u, v] + b^HF (k)[u, v] — e{u, v) , (76) 

with Win) and {{k) as before (but now on i/^/^(M^)), see ([701 . and 

b-yHF(K)[u, w] = a{u, Rj-hfv) — a{u, KjIif{k)v) , (77) 

where 

..,„w(x,.)=f ^^M-;7-7^ . (78, 

Using (i) of the proposition (exponential decay of the Hartree-Fock orbitals {^j}jLi) 
one now proves that (|78p extends to a holomorphic family of Hilberts-Schmidt op- 
erators in Sjy^ . One can now repeat the reasoning in the proof of (i) to obtain the 
stated exponential decay of ip. 

□ 

Liemma 5. (M, p. 196]) Let W{k) = e*''^ be a one-parameter unitary group (in 
particular, A is self- adjoint) and let D be a connected region in C with G -D. 
Suppose that a projection-valued analytic function P{k) is given on D with P(0) of 
finite rank and so that 

W{ko)P{k)W{ko)~^ = P{k + Ko) for kq G M and k, k + kq e D . 

Let ip G Ran(P(0)). Then the function ipii^) — W{k)^ has an analytic continuation 
from D r\R to D. 

Appendix A. Some useful lemmata 

Lemma 6. Let ft be an open subset of MP \ {0} with smooth boundary and let 
/i, /2 € H'^{ft) for some k >1. 
Then the function 

^ ' U |x-y| 

belongs to C'^(ri) if k>2, while if k — 1, it belongs to W^'P{fl) for all p > 1, and 
hence to C(f2). 
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Proof. We are going to prove the following equivalent statement. If A; > 2, xF S 
C'=(R3) for all X S C^(f^), while if fc = 1, G W^'P{M.^) for aU p > 1 and 

Fix X G C^{n) and take x G C^{Q) verifying x = 1 on suppx and such that 
there is a strictly positive distance between suppx and supp (1 — x)- We write 
xF(x) =xi^i(x)+ xJ^2(x) with 

^r(x)= / mflMIliyldy and F.(x) = / (1 - x(y))^^^f^ • 
Jr3 |x-y| Jr3 |x-y| 

The term x^2 is clearly in C°°(M"^). For the other term we use Young's inequality: 
if / e and g € L^R^) then 

||/*.9||. < with 1 + 1 = i + i. (79) 

r p q 

Moreover, ifl/p+l/g=l then f * g is continuous (see [24) Lemma 2.1]). Let 
a£Nl with \a\ < k. Then 

|i?"(xFi)(x)| < J2 \D'''xi^)\ I ^^c''nx/l/2)(y)dy . (80) 

If /i,/2 G -ff'=(f^), k>2, then i:'*(x/i/2) G L5/3(m3) for all /Ja as in From 
([75)1 . (|5D)) and x/l • I e L5/2(]^3-) follows that D"(xi^i) is continuous and, since a 
is arbitrary, that X-F G C'=(R3). 

If /i,/2 G H^{n) then 9(x,/i,/2) G ^^/^(M^) and from ([791) we get (only) that 
d{xF) e LP(M3) for all p > 1. It then follows that F e W^'P{n) for all p > 1 and 
therefore (by the Sobolev imbedding theorem) F G C{Vl). □ 

Lemma 7. Lei, for Za < 2/7r, /iq be the self-adjoint operator defined in (9), and 
let A_ (a) be the projection onto the pure point spectrum of Hq . 
Then the operator A-{a)hQA^ (a) is Hilbert- Schmidt. 

Proof. Let e > be such that Za{l + e) < 2/7r(l — e). We are going to prove that 
there exists a constant M = A/(e) such that 

with C Za{M + 2a-i)(l + 1/e) and P = X[o,m] (7'(p))- The claim wiU then 
follow from (|5T|l since 

^-([^"]-)^^ (M+L-i)^ ^-([-^-n^-)^<-- 

The last inequality follows since the eigenvalues of — A — C/ 1 • | are — / , n e N, 
with multiblicity n^. 

We now prove ([8T|l . For e > and any projection P (with P^ = 1 — P), we have 
that 

ho = PhoP + P^hoP^ ~ P^^P^ P^^^P 

>Piho-l^^)P + P^iho~e^^)P\ (82) 

By a direct computation one sees that there exists a constant AI — M{e) such that 
r(p) > M implies r(p) > (1 - e)|p| and T(p) < M implies T(p) > (-A). 
Hence, with this choice of M and P = X[o,a/] (^^(p)); l|82p implies that 
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The inequality ((ST|) follows directly by the choice of e. □ 

Appendix B. Pseudodifferential operators 

In this appendix we collect facts needed from the calculus of pseudodifferential 
operators ('0do's) (for references, see e.g. [5] or [H]). 

Define the standard (Hormander) symbol class S'^(]R"), /Lt G M, to be the set of 
functions a € C°°(M^' x M^) satisfying 

|5^5fa(a;,C)| < Ca,;3(l + l^n^^-l^'D/^ for all (x, e K" x . (83) 

Here, a, f3 e N" and \a\ = ai -\ h a„. Furthermore, 5^(M") C S'^'(M") for 

/i < Ai'. We denote S'°°(E") = U^6r5''^(R") and S'-°°(]R") = n^^RS''{W'). Finally, 
note that ab £ S'''i+''-^ (M"), 9^9|a e S'''!"!''! (M") when a G S''^i(M"),6 g ^'^^(M"). 

A symbol a e S'^(R") defines a linear operator A — Op(a) G: ('pseudodiffer- 
ential operator of order /i') by 

[Op(a)u](a;) = (27r)-" / e-«a(a;, 0^(0 , (84) 

where u is the Fourier-transform of u. The operator A is well-defined on the space 
iS(M") of Schwartz-functions; it extends by duality to 5'(R"), the space of tempered 
distributions. Note that for 

0<|a|<M 

(with smooth and with all derivatives bounded, i.e., Ua G B{M.'^)), A ~ Op(a) G 
^f** is the partial differential operator given by 

[Op{a)u]ix)= aa{x)D°'u{x) . (86) 

0<|q|<m 

Note also that, with a — a{x) and h — 5(^), 

[Op{a)u]{x) = a{x)u{x) and [d^u](0 = 6(0"(0 • 

If a e ^^(E"), then Op(a), defined this way, maps H^{W^) continuously into 
H^-t^i^^^) for all A; e M. Here, i/'^(M") is the Sobolev-space of order fc, consisting 
of M G 5'(R") for which 

ll"ll?^'=(E") := / \m\\^ + \etdi (87) 
is finite; this defines the norm on H^{W"). We denote 

iJ°°(R") = Pi i7'^(R") , i7-°°(R") = IJ i/'=(R") . 

feGK feGK 

In particular, symbols in S'°(IR") define bounded operators on L^(IR") = iJ''(M"). 
Furthermore, operators defined by symbols in S~°"{W^) maps any H''{W") into 
iJ°°(K"); such operators are called 'smoothing'. 

We need to compose ipdo's. There exists a composition ^ of symbols, 

# : S^"^ (R") X 5''"^ (R") ^ S^'i+A's (j^n-j ^gg^ 
(a,6)K^a#6, (89) 
such that 0p(a)0p(5) = Op{a^b). It is given by 

(a#5)(x,0 = -i- / e-*^-«a(x,^- 77)6(0: -y,r?)dydr?. (90) 

Here, the integral is to be understood as an oscillating integral. 
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The symbol a#6 has the expansion 



(91) 



a 



Here, 



means that for all j e N, 




(92) 



\o:\<j 



(recall that {d"a){d^b) e 5'*^+''^ l"l). One easily sees that the composition is 
associative. 

Proposition 2. If a e 5™i(M"), b e S'^'{W) then the symbol associated to 

[0p(a),0p(6)] belongs to S"''+"''-\W). 

In particular, if (j)i,4'2 G B°^{W^) (the smooth functions with bounded deriva- 
tives) with supp (^1 n supp ^2 = and a G 5''(M"), a{x,^) = a(^), then ^i:^a#(^2 ^ 
0, and so, with A := Op(a), 



is smoothing. 
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